Analytical results for the Coqblin-Schrieffer model 
with generalized magnetic fields 
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Using the approach alternative to the traditional Thermodynamic Bethe Ansatz, we derive an- 
alytical expressions for the free energy of Coqblin-Schrieffer model with arbitrary magnetic and 
crystal fields. In Appendix we discuss two concrete examples including the field generated crossover 
from the SU (4) to the SU (2) symmetry in the SU (4)-symmetric model. 
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The advances of nanotechnology have given an addi- 
tional weight to the old problem of Kondo effect. Ar- 
tificially manufactured structures such as quantum dots 
emulate the behavior of "natural" magnetic impurities 
though on different energy scales. Using technologi- 
cal means one can widely vary parameters of the dots 
thus getting an access to previously experimentally un- 
explored regions of the phase diagram. 

The most physically transparent situation corresponds 
to the case when magnetic impurity (or quantum dot) has 
a perfect symmetry. This, however, is rarely achieved in 
real systems due to a presence of the crystalline lattice. 
Let us consider, for instance, magnetic impurities made 
of rare-earth magnetic ions of Ce and Yb. In the presence 
of a strong spin-orbital coupling an f 1 (Ce) or / 13 (Yb)- 
orbital is characterized by the total angular momentum 
i (i = 5/2 for Ce and 7/2 for Yb) such that an isolated 
ion has the SU(n) symmetry with n = 2j + 1. In the 
crystalline environment this symmetry is broken. The 
interplay of these effects with the Kondo screening can 
be studied using the Coqblin-Schrieffer model 0: 
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where ct , c fc are creation (annihilation) operators of 
the conduction electrons partial harmonics with the an- 
gular momentum projection m = j + 1 — a (a = 1, . . . n); 
ft and f a -operators describe the impurity spin and V is 
the volume of the system. The generalized magnetic field 
H a originates from crystal fields inherent to the material 
and the external magnetic field. Since Hamiltonian l|T)l 
commutes with the operator 

n 
o=l 

one can assume without loss of generality that 

n 

E^ = °- ( 2 ) 



Notice that the cases described above correspond to the 
sector of the Coqblin-Schrieffer model with the occupa- 
tion number q = 1. 

Since model Q has a linear spectrum, it has to be 
equipped with the ultraviolet (UV) cutoff A and a con- 
sistent removal of the UV divergences requires that the 
"bare" coupling g a — njp(0) (here p(0) is the conduc- 
tion electron density of states at the chemical potential) 
be given certain dependence of the cutoff momenta (see, 
e.g. 0): 
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Eq.||3J) shows that for positive go the Coqblin-Schrieffer 
model acquires a physical energy scale, the Kondo tem- 
perature 



K 
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and renormalization trades the bare coupling constant go 
for the renormalization group invariant scale Tr. There- 
fore the partition functions of the model in the sector 
with a given occupation number q, Z q , actually depends 
on the dimensionless combinations T/Tk and H a /TK- Of 
course, formula J3J does not specify the physical energy 
scale uniquely and in order to define Tk unambiguously 
we shall impose the conventional normalization condition 
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where C\(T) = T-^(T log Z q= i) is the heat capacity of 
the impurity in the sector with the occupation number 
q = l. 

If the fields are weak in comparison with the Kondo 
temperature JSJ, the behavior of the system is governed 
by (the Kondo temperature of the fully SU(n)- 

symmetric model) which for rare earth impurities may 
be as large as several hundred degrees. On the other 
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extreme, if the generalized magnetic fields exceed 
and break down the degeneracy to one Kramers dou- 
blet, one gets the Kondo temperature of the order of 
several degrees. As an illustration let us estimate the 
new Kondo temperature for the case when the original 
SU(n) symmetry is broken down to SU(m) by the fields 

(n) 

^ Tpr . The estimate is easy if all fields 



K 

are of the same order H. Since InT 



(n) 
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— l/njp(0), lnT^™- 1 ~ — l/mjp(0) and the dimension- 
less ground state energy depends only on H a /T^\ we 
conclude that the resulting Kondo temperature is 
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which may easily constitute a scale vastly different from 
the Kondo temperature of the unperturbed model. It is 
clear that the detailed behavior in this crossover interval 
depends on the field ratios H a /Ht and it would be highly 
desirable to have analytic tools to handle a situation with 
an arbitrary pattern of fields. 

Historically the thermodynamics of the Coqblin- 
Schrieffer model has been examined by the method of 
Thermodynamic Bethe Ansatz (TBA) [1SHE0- Un- 
fortunately the TBA equations corresponding to Eq.QJ 
are rather complicated to be studied analytically, and 
until now the majority of results have been obtained by 
means of their numerical integration (see || and refer- 
ences therein). 

This work is based on the approach alternative to the 
TBA. Here we just outline major steps of our analysis 
and refer the reader to the papers |3, [lfj, [ill LL2J where 
the method itself was developed. Let us make some re- 
marks on the general spirit of the approach. It uses the 
fact that integrable impurity models in general and the 
Coqblin-Schricffer model in particular can be mapped to 
the (l-fl)-dimcnsional bulk conformal field theory (CFT) 
with a non-conformal boundary interaction (see e.g. |l3|). 
From this point of view, Eq.Q) belongs to the class of 
exactly solvable boundary theories such that the associ- 
ated boundary state commutes with the infinite set of 
mutually commutative local integrals of motion of the 
bulk system. As a matter of fact, one may say that the 
boundary state "generates" this set in the sense that it 
admits the asymptotic large-distance expansion in terms 
of these local integrals. More precisely, the Hilbert space 
of the bulk CFT associated with Eq.Q can be classified 
in accordance with the WA n _i- algebra with the cen- 
tral charge c = n — 1 and the corresponding boundary 
state commutes with the set of local integrals of motion 
introduced in [14| . 

In the approach adopted in this paper the amplitudes 
of the boundary state are related to monodromy charac- 
teristics of certain ordinary linear differential equations. 
The key ingredient is the equation of the form, 

{(- id v + hi) . . .(-id v + h n ) - e ne e"z;}* = , (7) 



where 9 and h a are some (complex) parameters. For 
3to(#) — 7r/2 the equation admits a solution which 
specified unambiguously by the following asymptotic as 
v — > +oo: 



( 
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Here C is an arbitrary constant, whose explicit value is 
not essential for our purposes. Notice that Eq.{7|) is in- 
variant under the transformation, 6 — ► 9+2Tri/n, whereas 
the asymptotic (JBJ is not. Hence the analytic continua- 
tions of ^0, 



%(v,6) = 9 (v,6 + 2iri q/n) , 



(9) 



with integers q, generate new solutions of Eq.fJTJ). 
It is possible to show that the Wronskian 
W[^0j •••^n-i] does not vanish, so that the 
set g is a fundamental system of solutions of 

Eq.0. By virtue of this fact, solution @ with q = n 
can be decomposed as 

n-l 

9) = ^(-l)"-"- 1 Z q {9 + inq/n) * q (v, 9) . (10) 

q=0 

Following the line of arguments similar to that of 0, El , 
it is possible to show that if the parameters 9 and h a 
are identified with the dimcnsionless parameters of the 
Coqblin-Schricffer model, 
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then the function Z q {9) appearing in (|lCj|> coincides with 
the analytic continuation of the partition function of this 
model, Z q , for the sector with the occupation number q. 
The subject of our current interest is the free energies 



T q = -T Iog(Z 9 ) 
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at the low temperature limit. In particular we study 
vacuum energies £ q — T q \T=o- They can be obtained 
from the semi-classical (WKB) approximation for Eq.Q. 
The leading terms in the WKB expansion of the solution 
*o © read, 



exp <^ 1 e 
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du (X(Y) - Y) 
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Here X = X(Y) is the solution of the algebraic equation 
Y n = {X + hx)...(X + h n ) , (14) 
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such that 



X(Y) -> Y as Y oo . 



(15) 



In fact, X = X(Y) is a multi- valued function of the ar- 
gument Y and Ea. H15(l uniquely specifies its branch for 
a sufficiently large Y. The latter solution admits a con- 
vergent 1/Y-power series expansion found by Lagrange 

M 



oo 
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Here Ik are symmetric polynomials given by 

(_]\c«i+...+a!„ 



Ik(hi, ■ ■ .hn) = 

r(ai + . . . + a„_i 
x i 

r(i - £) 
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G*2 1 G3 2 . . . G"™- 1 , (17) 



where we use the following notation for the elementary 
symmetric polynomials G&: 



l<ai<a2...<ak<.n 



(18) 



Note that according to the constraint (J2J, we have set 
G 1 = 0, and also G 2 = - ± E"=i h l- 

Using the Lagrange formula (| 1 fcil) and 113( 1 it is possible 
to show that the vacuum energies admit the following 
weak-field expansion |l6j 



£ q = Eo sin y — ) -I x 
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where 



G fc = 2 k+1 k k/n T k (l/n)T{-k/n) . 



The first term in l|19|) is the vacuum energy for zero fields 
and the value of En~ Tk is related to the choice of the 
constant C in (J""} 17\. Notice that Ik (|17J) are homoge- 
neous symmetric polynomials of the variables h a of the 
degree k + 1, i.e. 

I k (Xhl, . . . Xh n ) = X k+1 Ik(h!,... h n ) , (20) 

so Ea. l[T"""j> can be considered as a power series expan- 



sion in H/(2ttT k ), where H = ^E"=i ff «- Tllus if 
the fields are weak in comparison with the Kondo tem- 
perature the behavior of the system is governed by the 
convergent series ((19(1 . At the same time this series has 
a finite convergence radius and defines the multi-valued 



function of complex variable H (see Appendix for exam- 
ples) . 

We have also examined the structure of the low tem- 
perature expansion of the free energy ((12(1 and found the 
following asymptotic series: 



^ = W^W^ £ G fc sin(^) 
Hi H, : 
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Here the numerical coefficients Ck are the same as in 
(|19f) and Ife(/H, • ■ ■ , h n ) is a symmetric polynomial of the 
degree k + 1 of the form, 



h + Jk 



(22) 



where Ik is given by Eg. 1(17(1 and Jk = Jk(hi, ■ ■ ■ , h n ) are 
some symmetric polynomials of the degree k. Currently 
the polynomial tk are not known in a closed form for 
arbitrary k. Nevertheless they admit a simple algebraic 
description. As has been mentioned before, the bound- 
ary state associated with the Coqblin-Schrieffer model 
commutes with the infinite set of mutually commutative 
integrals of motion for the M / A„_i-algebra with the cen- 
tral charge c = n — 1. It turns out that the polynomials 
Ik appearing in Eq. I|21|) coincide with vacuum eigenvalues 
of these conserved charges. Therefore they can be calcu- 
lated purely algebraically for any finite k. In particular 
we found, 
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24 
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Here Gk are the elementary symmetric polynomials 11811 . 
Notice that for the physically interesting case of the mag- 
netic field configuration 



h a = h (j + 1 - a) , j 



n- 1 



(24) 



all polynomials Ik with even k vanish, whereas using 
Ed. 1(2*3} one has: 



h = ^ (n(n + l)h 2 + l), 

(n-l)(rc-3) 
1920n 

,2/ , ,W„ , n\ u4 



(25) 



(n 2 (n + l)(n + 3) h 4 + 10 n.(n + 1) /i 2 + 9 ) 



In conclusion let us just repeat the basic results of 
the paper. We have outlined a method which allows 
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to extract analytical results for the thermodynamics of 
Coqblin-Schrieffer model in the presence of magnetic and 
crystal fields. The general formulae drastically simplify 
for the SU(4) model (see Appendix). 
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Appendix 

In this appendix we illustrate the expansion (|19J) using 
two examples. 

Our first example is related to the case n = 4. Let 
us consider the following configuration of the generalized 
magnetic fields: 



Hi — —Hi , Hi 



-Ho 



(26) 



with Hi > H2 > 0. In this case Eo.fi 4jl is especially 
simple and can be written in the form 

Y 4 = (X 2 - h 2 )(X 2 - h 2 ) . (27) 

Its solution X = X(Y) satisfying 115|) reads explicitly: 



x = y \J yjl + t/y 1 + l/y 2 
where y = Y/h, x — X/h, and 



(28) 



h = 



h\ + h\ 



\h\ + h 2 ) 



All even coefficients I21 in series (|16f) vanish now, whereas 
I21-1 can be expressed in terms of the hypergeometric 
function: 



V ' sfiiW \ 2' 2 2' 



3_l_ 

4~2 ; 



Combining this equation with (|19|l one obtains, 
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with 
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\H 2 + H 2 ) ~ 



Eq. (I30|) can be rewritten in the form of a convergent inte- 
gral which is useful for an analytical continuation of the 
power series expansion outside its convergence disk, 



\H/T K \ < 



7TC 



V2F(l/4) 
Explicitly one has (q = 1,2,3): 
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It is particularly illuminating to extract from this ex- 
pression the Kondo temperature for the SU (2) Coqblin- 
Schrieffer model to compare it with Eq.©. The Kondo 

(2) 

temperature T K can be extracted from the magnetic 
susceptibility in the limit when one of the fields is very 
large. Since the energy depends on two fields H\, H 2 , 
one has to be careful in choosing the right direction of 
differentiation. The right choice of variables is H± = 
Hi ± H 2 such that d£i/dH„ = at = 0. Then at 
H+ ->■ 2H > T K , H- -> we obtain for q = 1 the 
following expression for the magnetic susceptibility 

d 2 £i 1 
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Hi — 2 H 
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\/8r 2 (i/4) .ff 



(32) 



The above expression for the effective Kondo tempera- 
ture is a particular case of Eq.© for n = 4, m = 2. 

Our second example is related to the case of an arbi- 
trary integer n. In |3j, the vacuum energy £1 was found 
for the following field configuration, 

'n(2a- 1) 



H a = V2 H cos (- 



2n 



(33) 



For this pattern the equation (|14|) is expressed in terms 
of the Chebyshev polynomials T n (x) = cos(n arccos(a;)): 



where y = Y/h, x — 
x= y 2F1 ( 

7T- 2-Fif— , 
2v V 2n 



y n = 2 1 -^ T n (x/V2), 
X/h. Therefore 

1 n- 1 



(34) 
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n 2n n 



,2n 



(35) 



and non-vanishing coefficients in (|16(l and (|19fl are 
/atJ-X (1 = 1,2...) and hnl+i (I - 0, 1, 2, . . .): 



Hrd+a 



r(f + 2i) 



nr(i 
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(36) 
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As in the first example, the vacuum energies l|19(l can be 
written in the form of convergent integral 



£ q = E sin ( — 
V n 

V2Hsin(TTq/n) f 00 ±c 



47TO 

V 2n 2 J 



2tt 



(37) 



T (±_^\ / nir(l/n)H 
\2n 2 M ny/2Tx 



It is easy to check that this equation for n — 4 is in 
agreement with (|3*l"1) provided e = 1/2. 
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